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-AEROX- 


COMPUTER PROGRAM FOR TRANSONIC AIRCRAFT AERODYNAMICS TO HIGH ANGLES OF ATTACK 
VOLUME I.- AERODYNAMIC METHODS AND PROGRAM USERS' GUIDE 

John A. Axelson 
Ames Research Center 

Summary 

The theory formulated in the AEROX computer program for estimating aircraft 
aerodynamics to high angles of attack (60°) and a program user's guide are 
included in the present Volume I. Program operators concerned specifically 
with program execution should consult pages 15 through 19. Volume II contains 
comparisons of the estimated and experimental aerodynamics for nine test cases 
over wide ranges of angle of attack and Mach number. A program listing and 
sample output tables and plots are shown in Volume III. 

The AEROX program estimates the coefficients of lift, induced-drag and 
pitching moment for wings and for wing-body combinations with or without an 
aft horizontal tail. Both trimmed and untrimmed characteristics are estimated. 
Minimum drag coefficients (e.g., friction, wave and propulsion-system additive 
drags) are not estimated in AEROX, but may be input as initialized coefficients 
for Inclusion in the total aerodynamic parameters, such as lift/drag ratio. 
Low-speed viscous stall is not covered. Transonic airfoil calculations utilize 
the chordwise locations of the local limit shock waves. These designated shock 
locations are entered as an input parameter, rather than extracted as solutions. 

The AEROX program is based on new forms of compressible wing theory 
covering potential and nonpotential flows, allowing for local regions of tran- 
sonic limit Mach numbers, and encompassing attached and detached leading-edge 
shocks in supersonic flow. The program is applicable to broad ranges of con- 
figurations and flight conditions, and has a rapid computer execution time, 

(30 points per second on the IBM 360). It has been used extensively in 
computerized aircraft preliminary design and optimization studies, and would 
be a valuable asset in activities concerned with aerodynamic instruction, 
correlation and research. 



INTRODUCTION 


The AEROX computer program for estimating aircraft aerodynamics to high 
angles of attack is based on new, nonlinear for lulations of compressible v/ing 
theory extended to cover transonic and supersonic flows, including shock waves 
and attendant rotational flows. The unavailability of any rigorous method for 
calculating three-dimensional aerodynamics to high angles of attack was dis- 
cussed by A.M.O. Smith in the 1974 AIAA Wright Brothers Lecture reported in 
reference 1. The primary obstacle cited was the continuing and current 
inability to define the overall influence of viscosity. 

The AEROX ap; roach postulates that transonic and supersonic flows past 
airfoils are dominated by shock waves, aud that the primary effects of viscosity 
are in determining the skin friction and in influencing the chordwise location 
of transonic shock waves. These shock locations are not extracted as solutions, 
but rather are entered as a designated input parameter. The effects of varying 
the shock location are then visible in the listed and plotted outputs of aero- 
dynamic coefficients. 

The nonpotential lift and induced drag equations are based on the inte- 
grations of momentum components in the plane through the trailing edge and 
normal to the flight direction. The onset of nonpotential, transonic flow 
around conventional, blunt airfoils is assumed to occur when the Laitone limit 
values for local Mach number or loading from references 2 or 3 are attained on 
the airfoil nose or along its surface. Additional discussion and supporting 
experimental evidence for this concept appear in reference 4. 

The prototype program has been in operation at the Ames Research Center 
since 1973, (FORTRAN IV, LEVEL G) , on the IBM 360 and CDC 7600 computers. It 
has been used in numerous aircraft preliminary-design and optimization studies, 
some examples of which are documented in references 5, 6 and 7. The synthesis 
program also included the CONMIN optimizer repaired in reference 8. 
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Nomenclature 


The symbols used in the input-output formats, the data set listings and 
in the equations are identified as follows: 


AEROX 

ALEL J, J 

ALER, CL 
h 


ALFTR * “prim 

ALPHA, ALF, a 

ALTV 

AMC 


APLAN 

ARDET 

ARH 

ARW 

ARWX 

ASECT 

AT62 

AXE 

BCLCD1 

BDMAX 

CBARW 

CDCAM 

CDHOR 

CDN 

CDO 

CDOB 


CDSEP 
CDTOT, CD 
CDW, CDj 
CISQ 

CLALFA 

CLHOR 


program name 

input integer identifying type of airfoil (14 J f- 5) See AXE listing 

angle of attack for onset of the leading-edge limit Mach number, 
rad. 

trimmed angle of attack, deg. (ITRIM=1) 

angle of attack of wing reference plane, deg. 
input altitude, ft. 

maximum angle of attack for subsonic compressibility or 
interference (40) deg. 

planform area of nose, sq. ft. 

angle of attack for shock detachment from sharp airfoils, rad. 

input aspect ratio of horizontal tail 

input aspect ratio of wing 

aspect ratio of exposed wing 

nose maximum cross sectional area, sq. ft. 

subroutine for atmospheric properties 

subroutine for wing aerodynamics 

subroutine for body aerodynamics 

input body diameter, ft. 

wing mean aerodynamic chord, ft. 

drag coefficient due to wing camber CLO 

induced drag coefficient of horizontal tail (ref. to SWING) 

nose or body induced drag coefficient (ref. to SWING) 

input minimum drag coefficient for wing (exclude camber) ICD0=1. 

input additive drag coefficient for body, tail, propulsion, 

I CDO 1. 

wing separation drag coefficient in flow zone 4 

total drag coefficient, includes CDO, CDOB when ICDO“l. 

wing induced drag coefficient 

polynomial in wing lift-curve slope equation, 
sinWcos^-zsin 2 **- £ SIN^COS 2 ^ 

wing lift-curve slope, per deg. 

horizontal tail lift coefficient (ref. to SWING) 
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Nomenclature - Page 2 



CL1.ES 

CLN 

CLO 

CLOB 

CLTOT , CL 

CLU 

CLWL 

CLWU 

CLUBS 

CLUES 
CM, C 

m 

CMO 

CMOB 

CPDS 

CPLIM 

CPMXS 

CPSTAG 

CP 2 

CROOT 

CTIP 

DALTR 

DCDTR 

DCLTR 

DDISP 

DELH 

DEPDA 

DEXP 


DLWING ~ 

’dc< 


DNBL 

DWASH, EPS1 


wing lower-surface lift coefficient at a 

£ 

nose or body lift coefficient (ref. to SWING) 
input wing lift coefficient at a=0 

input additive lift coefficient for body, propulsion 
total lift coefficient 
wing lift coefficient 

lift coefficient fox wing lower surface 
lift coefficient for wing upper surface 

incremental wing lift coefficient of upper surface in interval 
aH *E 

wing upper surface lift coefficient at a 

£ 

pitching-moment coefficient 

input wing pitching-moment coefficient at a=0. 

input additive pitching-moment coefficient for body (BCLCD1 
subroutine) 

pressure coefficient downstream of limit shock for Z=4 (IXCD=0) 
limit pressure coefficient 

pitot compressibility factor for lower surface lift (Zg4) 
stagnation pressure coefficient behind normal shock (Z=6) 
pressure coefficient downstream of limit shock for Z=4 (IXCD=1) 
wing root chord ( 4 ) , ft. (streamwise) 
wing tip chord, ft. (streamwise) 

adjustment of angle of attack for trim at constant C , deg. 
incremental trim drag coefficient 
increment of lift coefficient during trim 

plot control integer for component drag coefficients; 0, No plot, 
1, plot printed 

increment of horizontal-tail deflection to trim, deg. 

rate of change of downwash angle with angle of attack 

exponent in induced drag coef f icient equation, 1.5-Mcos SWPWLF 
for M 1; 1.5-M 2 cos SWPWLF for 1 M. * 

4 

wing lift-curve slope, per rad. 

Mach number attenuation factor for lift-curve slope in Z*6 
downwash angle at the designated horizontal tail location, deg. 
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Nomenclature - 
ED 

EPSD 

FDNOSE 

FINT 

FLAX 

FLEX 

FLNOSE 

FMOML 

FMOMU 

FTOTD 

Fl'OTL 

FCM 

GEOM1 

ICDO 

IDATA 


Page 3 

eccentricity of wing CBARW/4; XQMAC-XCG 
wake downwash angle 

multiplying factor for nose drag, programmed unity 

wing-body interference factor, attenuated with a, Z<4, M<1.4 

wing-body interference facto-, 5<Z 

lift factor for wing chord extension 

multiplying factor for nose lift, programmed unity 

wing lower-surface lift factor for Z=6 

wing upper-surface lift factor for Z=6 

multiplying factor for change of reference area for 

multiplying factor for change of reference area for C 

L 

multiplying faccor for change of reference dimensions for C w 

M 

subroutine for calculating geometric parameters 

input control integer for minimum drag; 0, CDO omitted; 1, input 
wing CDO included 

input control integer on printout; 0, no print; 1, listing 
printed 


IFLEX 

IPLOT 

1TABL 

ITRIM 

IT 

IXCD 

JjALEU 

L/D 

LDISP 

LE 

LT 

M, SMN 
Me 

MPLOT 

P 2 

PLOT 

PPLOT 


input control integer for strake bluntness; 0, sharp; 1, blunt 

input control integer for plots; 0, no plots; 1, plots printed 

input control integer for geometry heading; 0, no print; 1, 
heading printed 

input control integer for trim option; 0, data untrimmed; 1, 
trimmed 

input horizontal-tail incidence, deg. 

input control integer for limit shock position; 0, constant X/C; 
1, limit shock sweep angle SHK specified from XCD at airplane 
centerline 

input integer identifying type of airfoil - See AXE listing 
lift/drag ratio, when ICDO-1. 

input plot-control integer; 0, no plot; 1, component lift co- 
efficient plotted 

tail length from moment center, ft.; XQHOR-XCG 
tail length from CBARW/4, ft.; XQHOR-XQMAC 
Mach number 

cross-flow Mach number, M sina in BCLCD1 subroutine; fig. 3. 
subroutine for trimmed plots, all Mach numbers together, 
static pressure downstream of limit shock 
subroutine for untrimmed plots, plot set for each M. 
input plot-control integer; 0, + Cm right; 1, + Cm left 
free-stream stagnation pressure 
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Nomenclature - Page 4 


RNC 

RHLOC 

ROC 

SEXT 

SHK 

SHOR 

SMN,M 

SPANW 

SQH 

sqwjsvri'w 

SWING 

SWPWLE 

SX 

SXCjSHK 

TCRW 

TCTW 

TRIMDG 

TRW 

XCD, X/C 

XCG 

XEXT 

XLB 

XLN 

XQHOR 

XQMAC 

Y 

YBAR 

YHOR 

Z 

a, ALPHA 
dg, ALER 

Y 


body cross-flow Reynolds number; fig. 3. 

Reynolds number per foot for input Mach number and altitude 
input leading-edge radius to chord ratio for J=5 airfoils, 
input area of wing chord extension (forward strake), sq. ft. 
input sweep angle of the limit shock in Z=4 with IXCD=1, deg. 
input horizontal tail area, sq. ft. 

Mach number 
wing span, ft. 

input sweep angle of horizontal-tail C/4 line, deg. 

input sweep angle of wing C/4 line, deg.-, rad. 

input wing reference area, sq. ft. 

Sweep angle of the wing leading edge, rad. 

exposed wing area, sq. ft. 

sweep angle of limit shock, Z=4, 1XCD=0;1. 

input thickness-to-chord ratio of wing <£. (streamwise) 

input thickness-to-chord ratio of wing tip chord (streamwise) 

subroutine for downwash, tail aerodynamics and trim 

input wing taper ratio, CT1P/CR00T 

input chordwise location of limit shock, Z=4. 

input longitudinal station of moment center (or CG), ft. 

input longitudinal station of centroid o f wing chord extension 
SEXT, ft. 

input body length, ft. 
input nose length, ft. 

input longitudinal station of horizontal-tail C/4, ft. 
input longitudinal station of wing CBARW/4, ft. 
vertical coordinate 

lateral station of the wing mean aerodynamic chord, ft. 

input height of horizontal tail from wing plane, positive for 
high tail, ft. 

integer identifying flow zone 

angle of attack of wing reference plane, deg. 

angle of attack for onset of the leading-edge limit Mach number, 
rad. 

ratio of specific heats for air, 1.4. 
horizontal-tail angle of attack 

bar over a parameter denotes the effective value; e.g., 

SQW - SIN” 1 [SIN(SQW)COS(a)] 
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AERODYNAMIC THEORY 


The equations for evaluating airfoil aerodynamics are organized according 
to the division of the flight envelope into the flow zones depicted in figure 1. 
The cross-hatched viscous stall region is not included in the present method. 

The flow is considered Incompressible for Mach numbers below 0.1. The com- 
pressible, shocklesu zone 2 covers airfoils with L'.unt leading edges and with 
surrounding flows having local Mach numbers everywhere below the Laitone limit 
value (SQRT(Y+3)/2) . The boundary between zones 2 and 3 is estimated by the 
equation for ALER, the angle of attack for onset of the limit Mach number at 
the leading edge. The onset of the surface limit-Mach number, zone 4, occurs 
when the upper-surface lift for zone 2 or 3 reaches the limit lift corresponding 
to the designated chordwise shock location, XCD. Sharp airfoils are considered 
to have reached the leading-edge Mach limit for all subsonic flow, and are 
treated with the nonrotentlal lift equations of zones 3 and 4. Zone 5 applies 
only to sharp airfoils having supersonic leading edges with attached shocks. 
Finally, flow zone 6 covers all elrfoils with detached, leading-edge shocks. 

The supersonic leading edges for zones 5 and 6 occur when the normal component 
of Mach number exceeds unity (M cos SWPWLE^l . ) . 

The following equations apply to an equivalent wing having straight 
leading and trailing edges. Account is made for strakes and wing forward- 
chord extensions in AEROX by multiplying all of the following equations for 
lift coefficient and wing lift-curve slops by the term (1 + FLEX). FLEX, the 
empirical lift factor for the chord extension, is defined within the program 
and depends on Mach number, angle of attack, and whether the extension is 
sharp (IFLEX-0) or blunt (IFLEX-1). 


Incompressible Flow Zone 1 (M^O.l) 


The incompressible lift equation for conventional airfoils is taken from 
the potential-flow theory of Kutta-Joukowski. The lift equation for nonpotential 
flow, such as around sharp airfoils, is based upon the integration of downwash 
momentum. Both equations are extended to three-dimensional flows through the 
inclusions of the Prandtl aspect-ratio transformation and the cosine term 
involving the effective sweep of the quarter-chord line. The potential-flow 
equations for conventional airfoils are: 

C L * 27f sinoc cos SQW (^j) (1) 


dC L 

j— - = 21T cosot cos SQW 
a \ 


c 


ARW 

ARW+2 


) 




IT ARW 


( 2 ) 

(3) 


The nonpotential equations derived in the Appendix are: 
C L - 2* S i,«o S 2 «(l-2£»)„ 5 sw( 7 ^) 


(4) 


dC L 4 2 2 2 2 ARW 

n — =1Tcos« (2sin ot+2cos «-4sin <x-3sin otcos <x)cos SQW ( -- -■ _) (5a) 

dot ARW+Z 

- 2ircos*(CISQ)cos SQW (j—j) (5b) 


C Dj - C l (TAN ) DEXP (6) 

The incompressible lift is comprised of equal contributions from the upper 
and lower wing surfaces. The nonpotentinl equations asymptotically approach 
the potential equations at small angles of attack. 


Compressible, Shock-free Flow Zone 2 

In estimating the aerodynamics for airfoils having subsonic leading edges, 
separate compressibility factors are ..ced for the lifts of the upper and lower 
surfaces. For compressible, shock-fre* flow (zone 2) on the upper surface of 
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blunt airfoils, a Prandtl-Glauert factor is used involving the component of 
flight Mach number normal to the quarter-chord line and an angle-of-attack 
attenuation to value unity at 40°. 

z 


(i-Vi-m 2 co 3* 5qvy)fjgj 
yi-M^cos 2 SOW 


No compressibility factor is used for the uppe^-surface lift in nonpotent. al 
flow, because the onset of the local Mach-li T .it is considered to "freeze" the 
local flow. 


Lower-surface lifts in potential and nonpotential flows arcund airfoils 
having subsonic leading edges are evaluated using a pitot-compressibility factor 
expressing the ratio of the leading-edge, stagnation-line pressure in compres- 
sible flow to that in incompressible flow. Account is included for variations 
in the angles of attack and sweep, aud for wing-body interference patterned 
after the FLAX parameter from references 9 and 10. 


F, = CPMX5 - 


fM COS 5WPWLE 


J + (^l)(F, N T\M« S , SWPWie] -it (*) 





(M«f) 

(So) 

(kMtl.4) 

(8b) 


(»c) 


The compressible, potential-flow equations for conventional airfoils are 
obtained by incorporating the compressibility factors, equations (7) and (8), 
into equations (1) and (2). 
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( 1 ) 

( 10 ) 


M 


) 


C L* [fU f| NT) cos saw 


r - w l 

V Tm 


( 3 ) 


Sharp aiifoils in subson-'c, compressible flow are treated as nonpotential, 
Mach-limited flows, and assigned to t]ow zones 3 or 4. 

Leadirg-edR2 Mach-.limited Flow Zone 3 


The onset of the Laitone limit Mach number in the curvilinear flow around 

a blunt-airfoil leading edge is estimated by the value of the onset angle of 

attack, cx_, derived as equation (All), in the Appendix. A brief discussion of 
r* 

the conceptual flow model undergoing the transition from compressible, potential 
flow to MaJn-limited, nonpotential flow follows. 


After attaining the local limit Mach number around the leading edge, 
further increases in angle of attack or in flight Mach number tend to enlarge 
the local Mach-limited region on the airfoil nose. The enlargement of the 
constant velocity profile In the curvilinear i ow around the nose disrupts the 
radial equilibrium between the local, radial static-pressure gradient and the 
local, centrifugal forces. (Irrotational, curvilinear flow exists only when 
radial equilibrium prevails ,whicr. requires an essentially inverse relationship 
between local velocity and local streamline radius of curvature.) The constant- 
velocity profile may be accompanied by the dominance of the local centrifugal 
forces and the possible formation of a separation bubble. When the locally 
supersonic flow decelerates (and reattaches) just downstream of the nose (with 
or without the separation bubble), a "peaky" pressure distribution results, 
(shown in ref. 4). This flow is designated in AEF.OX as leading-edge Mach- 
limited, or zone 3. When the separation bubble la present, the mixing action 
of the rotational-flow layers emenating fiotn the disrupted nose flow may promote 
the flow reattachment. 
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When the supersonic limit-Mach number extends well back on the airfoil 
upper surface, resulting in the flat pressure distribution (also in ref. 4), 
the flow is surface Mach-limited, and assigned to zone 4. Surface Mach-limited 
flow is modeled to have separation downstream of the surface limit shock, 
resulting in an additive drag component, CDSEP. Experimental confirmation of 
the separation appears in pressure distributions, samples of which were 
included in reference 4. 

Sharp airfoils in compressible, subsonic flow are treated as leading-edge 
Mach-limited flow with the nonpotential zone 3 equations. 

Ci-H)*™ «&('- ‘ 4 *)“* *w(a&) 

&0 = Cl (tan ecJ 0eXP 


M 

M 

(») 


The lift ci rves for blunt airfoils traversing from compressible, potential flow 
(zone 2) into zone 3 are assigned a continuous mathematical transition, rather 
than a discontinuous jump, so that the program may be coupled to an optimizer 
program. The lift equation is the sum of the zone 2 lift (eq, (9)) evaluated 
at o(g and an incremental lift for the interval («-©£) in zone 3, until ec|. M is 

reached. 


C L * Greater af 



[F V (FIN1)+F,]*sin * E <»sSJW 

£*.( 11 ) 



Surface Mach-limited Zone 4 


The conceptual flow model depicting surface Mach-number limited flow past 
an airfoil is shown in figure 2. The Laitone limit Mach number (refs. 2, 3) 
corresponds to the maximization of the static pressure behind the surface limit 
shock. 





t 


This criterion is extended to swept wings in the AEROX program with the 
following limit pressure coefficient. 



The airfoil lift coefficient for zone 4 is the sum of the lower-surface 
nonpotential lift (eq. (11) and the upper-surface limit lift, which is 
dependent upon the designated chordwise location, XCD of the surface limit 
shock. The program uses XC = f(XCD) which moves the limit shock to the trailing 
edge as the sonic leading-edge condition is reached. 

C L a - [ Cf L»M Cp DS ^ “ xc )] c <>s Fj( irsinacosVfl-^Jcoss^ (mv^) i7 ^ 

5~ = [ c W xc ) + 2 c P05 cosoc (cisq)cos 0*) 

v- dsw +c ‘-( ton *) +Cd *p (») 


CPDS is the pressure coefficient downstream of the limit shock, corresponding 

to p 2 in figure 2. The separation drag coefficient, C , accounts for the 

°SEP 

momentum deficit in the modeled separation wake passing the trailing edge. 

The wake is assumed to originate at the base of the limit shock, to have its 
upper edge follow a line inclined at one-half the angle of attack, and to have 
a linear velocity profile between the zero value on the surface to the free- 
stream value at the upper edge. 


Cp D5 


z 

m* 


J 


J 

£r|i 

1 

i 


i 


f+i 


(v-l) (^?)+M4IW 1 5XC4“2 


r 

TT 


-I 






Supersonic Attached-Shock Zone 5 


Airfoils having sharp leading edges with attached shocks are treated in 
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zone 5. The upper boundary is the angle of attack for shock detachment, ARDET, 
which is evaluated by equations based on curve-fitting figure 4 in reference 11 
for supersonic speeds, and based on the theory of reference 12 for hypersonic 
speeds. The supersonic lift coefficients are estimated by the nonpotential 
equation (11) using the exposed wing area, an empirical, Mach-attenuated aspect- 
ratio transformation, and the wing-body interference factor of Flax (ref. 9, 10) 
applied to the lower surface lift. 

H*) (2C) 

M' 

2 - 1/2 

For multiply the above equations by (M -1) 

At Mach numbers above 3, lift coefficients are evaluated by the explicit, 
oblique-shock theory of reference 12. 



Cps C L tan* 
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Supersonic Detached-Shock Zone 6 


M 


! 


Blunt airfoils with supersonic' leading edges and sharp airfoils with 

detached leading-edge shocks are treated with the following equations combining 

the nonpotential equation (11) derived from diverted momentum and modified 

Newtonian impact theory, which assumes increased importance on the windward 

surface loading at hypersonic speed. The limit pressure coefficient for the 

upper surface undergoes a transition from a value of -0.68 at M=1 and approaches 

the value - —2, the well-known empirical value suggested by MAYER in reference 13. 
M 



M 2 ' M” 



The upper surface lift coefficient becomes a small fraction of the total lift 
coefficient as Mach number increases. 


-c, 


cos« 


C Lv - Lessee OF < 




ARWX 
ARWX + 


ft A 


(26 a, 6 ) 





Cl* c l u + Cljj . 


dCfc_ f**.*Lj 
do( d X ao< 


C Dj r C^Tdrux 



The empirical constant? entering the diverted-momentum (nonpotential) lift 
contribution include DNBL, establishing the variation with Mach number, and 
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FMOMU and FMCML, expressing the division of the diverted momentum lift to the 
upper and lower wing surfaces. 




DNBl =• < 




FMOMU = 

< , 


X 1. 


M-.Z 


M- 1.77) 
(l.77^l y 0 

Hilt) 

(l.Z< M) 


FMOML= 


I.6-.6M 

0.7 

1 * M- .55 




APPLICATION TO AIRCRAFT CONFIGURATIONS 


Input 


Input file formats for the IBM 360 and CDC 7600 computers are enclosed. 
Each item is discussed here and also identified in the nomenclature. Note 
that parametric studies may be performed by entering the minimum and maximum 
values and increments for the wing aspect ratio (ARW), wing taper ratio (TRW), 
wing quarter-chord sweep angle (SQW), angle of attack (ALF), and the designated 
chordwise shock locations (XCD) . For a particular airplane or single value 
to any of these parameters, enter the value as the minimum and the maximum, and 
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NJ 


) 


assign an appropriate non-zero value for the increments, (see sample input 
sheets for test cases in vol. II). The program allows an array of up to(l) 
twenty Mach numbers and up to (K) twenty angles of attack. Angles of attack to 
40° are usually input as 2° minimum, 40° maximum, with 2° increments. Angles 
to 80° can be input as 4° minimum, 80° maximum, with 4° increments. Maximum 
angles should not reach 90°. Enter dot^-set title f (sec input forrv\cctj. 

NSMN number of Mach numbers in the array, up to 20. 

SMN list the Mach numbers in the array. 

ICDO input control integer for minimum drag; 0, CDO omitted; 1, input 

value of CDO included in CDTOT and L/D. 

CDO input CDO values for each Mach number in the array. 

CMO input CMO values for each Mach number in the array. 

CLOB input CLOB values for each Mach numb e r in the array. 

CDOB input CD0B values for each Mach number in the array. 

CMOB input CMOB values for each Mach number in the array. 

ITRIM input control integer for trim option; 0, untrimmed; 1, trimmed. 

ALELJ input airfoil identification integer; 1, sharp; 2, NACA 230XX and 
00XX airfoils; 3, MCA 6 series airfoils; 4, slab airfoils with 
round leading edges; 5, leading-edge radius-to-chord ratio specified. 

MNARW minimum wing aspect ratio 

MXARW maximum wing aspect ratio 

INARW incremental wing aspect ratio 

MNTRN minimum wing taper ratio 

MXTRW maximum wing taper ratio 

INTRW incremental wing taper ratio 

MNSQW minimum wing quarter-chord sweep angle, deg. 

MXSQW maximum wing quarter-chord sweep angle, deg. 

INSQW incremental wing quarter-chord sweep angle, deg. 

SWING wing reference area for equivalent wing having straight leading and 
trailing edges, sq. ft. 

TCRW streamwlse thickness-to-chord ratio at wing centerline 

TCTW streamwlse thickness-to-chord ratio at wing tip 

ROC leading-edge radius-to-chord ratio; specify for ALELJ"5; otherwise, 

zero, 

SEXT area of wing chord forward extensions or strakes, sq. ft. 

XEXT longitudinal station of centroid of chord extensions, ft. 

IFLEX input 0 for sharp chord extensions, 1 for blunt. 
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BDMAX maximum body diameter, ft. 

XLN nose length, ft. (longit. 0 station at nose) 

XLB body length, ft. 

XCG longitudinal station of moment center or center of gravity, ft. 

SHOR horizontal tail area, sq. ft. 

XQHOR longitudinal station of horizontal tail C/4 line, ft. 

ARH aspect ratio of horizontal tail 

SQH sweep angle of C/4 line of horizontal tail, deg. 

YHOR vertical height of tail from wing plane, ft. (+ tail above plane) 

IT horizontal tail incidence, deg. 

MNALF minimum angle of attack, deg. (usually 2°) 

MXALF maximum angle of attack, deg. (usually 40°) 

1NALF incremental angle of attack, deg. (usually 2°) 

MNXCD most forward chord location of limit shock in Z=4, (typ. .3) 

MXXCD most rearward chord location of limit shock in Z^, (max 1.) 

INXCD incremental fraction of chord for intermediate shock locations. 

IXCD control integer; 0, limit shock in Z=4 is at constant fraction of chord 

across the wing span; 1, limit shock starts at XCD at , extends 
outboard at the specified limit shock sweep angle, SHK. 

SHK limit shock sweep angle, deg. for IXCD=1; 0 for IXCD=0. 

ALTV altitude, ft. 

FTOTL multiplying factor for change of lift coefficient reference area (def.l.) 

FTOTD multiplying factor for change of drag coefficient reference area (def.l.) 

FCM multiplying factor for change of pitching-moment coefficient reference 

dimensions 

IDATA listing-control integer; 0, no print; 1, printout provided. 

ITABL listing-control integer; 0, no heading; 1, heading printed. 

1PL0T plot-control integer; 0, no plot; 1, plots provided. 

PPL0T 0, + CM plotted to right; 1, + CM to left. 

LDISP component lift-coefficient plot control integer; 0, no plot; 1, plot 
provided. Six outputs are CLHOR, CLN, CLWL, CLWU, CLW, CLTOT 

UDiS? component drag-coefficient plot-control integer; 0, no plot; 1, plot 
provided. Four outputs are CDN, CDSEP, CDW, CDTOT 

In addition to sweep angle, SQW, at least three additional wing planform para- 
meters must be specified: SWING, ARW, and TRW. The program calculates the 
following parameters if they are input as zero: SPANW, CROOT, CTIP, CBARW, 

YBAR. Also, for a wing alone having longitudinal station 0, at the leading-edge 
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apex, the program calculates XQMAC, if it is input as zero. For all other cases, 
XQMAC must be input. The choice of horizontal tail area depends on the tail 
configuration. Use total planform area for large tails on small afterbodies. 

With wide bodies enclosing twin jets, use exposed tail area. 

AEROX includes the following subroutines: 

1 . AEROX 

2. AT62 

3. GE0M1 

4. AXE 

5. BCLCDl 

6. TRIMDG 

7. PLOT 

8. MPLOT 


Driver control, read and write statements 
Atmospheric properties 
Geometric parameters 

Wing aerodynamics and configuration totals 

Body aerodynamics 

Tail aerodynamics, C„, trim 

M 

Un trimmed vs. a, vs. Cp, vs. for each M 

Trimmed C, vs. a, vs. C^, all M's together 
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COLUMN Z 

l 

T 


ism 3 60 In put Format 


TITLE up TO 56 CHARACt1R$ LONG 
fc ARRAYS 
NSHN • 

■ 

icon • 
coo • 

- chd • ' ~ 

CLOB • 

COOB • 

~ "CWOB iT 
ITAIM • 

«ND_ 

GWINGIN 

Aieu • * 


NNARW • 


f 

MXARW • 

9 

Inarm ■ 

f 



NNtrW • 


f 

mxtrW ■ 

9 

~Intrw • 

9 



MNjOW • 


9 

MXSOM • 

9 

INSOM • 

9 

CUV • 

9 

S«!NG ■ 


9 

SPANN ■ 

9 

CROOT • 

9 

CTIP • 

9 

TCAW • 


9 

TCTW V 

9 

xw*e~* 

9 

CIARW. 

9 

ROC • 
SEND 


9 

SEXT • 

9 

XEXT • 

9 

IPLEX • 

• 

wosfn* 









bdmax « 
(End 


9 

XLN • 

9 

XLB « 

9 

XCO m 

9 

Stailjn 
SHOR ■ 


9 

XOMOR • 

9 

ARH • 

9 

SON - 

9 

YHOR • 


9 

IT ■ 

9 





“"SEWfiT 









gflowin 
mnalf ■ 


9 

mxalf • 

9 

INALP • 

9 



HWxCO • 


9 

NXXCO ■ 

9 

!NXCB _ ft 

9 



IXCD • 


9 

SHK • 

9 





altv • 


9 







kpnO 

GFXCTOR 









ftotl • 

SEND 

SOofPUT 


9 

FTOTD • 

9 

PCM • 

9 











IDatA • 


ITABL • • 






-' I plot • " 

V 

mot • t 






LDJSP • 

i 

9 

9 9 9 

9 





DOI sP • 
‘ SEND 

• 

9 

9 9 



















- 
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CDC_7c3 0 iliPUT FORMAT 

Columns 

i 

0000 1 00 TITLE up TO 56 CHARACTERS LONG 
0000200 iARRAYS 


0000300 

NSmN 

■ 

t 








0000400 

smn 

a 

9 

9 

• • 9 

9 

9 

( comma 

NUMBER 

20 

0000900 

JCDO 

a 

f 

9 

• • • 

9 

9 

w 

a 

a 

0000600 

COD 

a 

9 

9 

• • • 

9 

9 

M 

n 

a 

0000700 

CMO 

a 

9 

9 

Ml 

9 

9 

•I 

H 

a 

0000800 

CLOB 

a 

9 

9 

f 0 • 

9 

9 

M 

»* 

a 

0000900 

CDOB 

a 

9 

9 

• • • 

9 

9 

H 

a 

a 

000 1 000 

CMOB 

a 

9 

9 

t • • 

9 

9 

H 

a 

a 

0001)00 

ITRIH 

a 

9 

9 

• • • 

9 

9 

M 

a 

a 


000) ?00 SEND 
0001 300“ JWlNGIN 
0001400 ALglJ • 


0001500 

mnarw ■ 

9 

mxarw 

a 

9 

InARW 

a 

9 




0001(00 

MNTfttf • 

9 

HXTRW 

a 

9 

fNTRW 

a 

9 




0001700 

MN$QW ■ 

9 

MXSOW 

a 

9 

INSOW 

a 

9 

cm 

a 

9 

0001800 

SWING * 

9 

SPANW 

9 

9 

CROOT 

a 

9 

CTIP 

a 

9 

0001900 

TCrw • 

9 

TCTW 

a 

9 

XQMAC 

a 

9 

CBARW 

a 

9 

0002000 

ROC ■ 

9 

SgXT 

a 

9 

Xf XT 

a 

9 

iflex 

a 

9 

0002100 

*EMD 











0002200 

(NOSE IN 











0002300 

BDNAX * 

9 

XLN 

a 

9 

XLB 

a 

9 

XCG 

a 

9 

0062400 

(END 











“ 0002 500' 

(YAIIIN 











0002(00 

SHOP ■ 

9 

XOHOR 

a 

9 

ARH 

a 

9 

SOH 

a 

9 

0002700 

YHOR • 

9 

IT 

a 

9 







0002800 

send 











0002900 

(FtOWIN 











OOO 3 OOO 

nnalf • 

9 

HXALF 

a 

9 

I NAIF 

a 

9 




“0003100“ 

mn^cd ■ 

9 

MXXCD 

a 

9 

INXCD 

a ' " 

9 




0003200 

IXCO • 

9 

SHK 

a 

9 







0003300 

ALTV ■ 

9 










0003400 

send 






. — . 





OOO 35 OO 

(factor 











0003(00 

FTOTl • 

9 

FTOTD 

a 

9 

FCH 

a 

9 




“"0003 700 

“IlfiD 





- 



. _ 

~ r_ ’’ ’ 4 



OOO38OO (OUTPUT 

000)900 (DATA • * ITABL ■ • 

0004000 " I PLOT ■ 1 PPLOT • . 

0004100 LDISP • t t * f t • 

0004200 001 SP • • • . t » 

“ 0004S00 afno “ 

Nine odditiona] curds, siariirxy Q~t column 7 } ore required for CDC 7600 operation.: 

At front of deck : pp06RAM AERoX KJPU T? OUTPUT, TAPES = INPUT } TAPE 6 * OUTPUT) 

At Tndlfde cK:~ TuHCVOfTi iVRSiN (A) 

ARSIN » A Siw(A') 

RETURN 

ENO . . 

Function cotan (a) 

COTAkj = i./tan(A) 
return 

END 
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APPENDIX A - WING (AXE subroutine) 

Derivation of Nonpotential Lift Equation 

A description of the conceptual transonic-flow airfoil model (fig. 2) 
involving local limit Mach numbers, shock waves, rotational flow and separation 
is included in the text under flow - zone 3. The lift no longer conforms to 
potential theory, but rather is evaluated in AEROX by equations based on the 
integration of the downwash momentum in the plane passing through the airfoil 
trailing edge and normal to the flight direction. The distribution of the 
downwash velocity in the vertical plane is defined by the nondimens ional 
attenuation factor appearing in the brackets of equation (A2) below. This 
factor prescribes the distribution of the disturbance velocity in the XT plane 
for rectilinear flow with superposed circulation (vortex). It also closely 
represents the distribution of the downwash velocity in the same XT plane also 
passing through the trailing edge of the cambered, streamline airfoil shown in 
the upper sketch. This factor also prescribes the distribution of the distur- 
bance velocities in the two components in the XT plane through the trailing 
edge of the planar airfoil shown in the lower sketch. The trailing edge 
velocity Vcosa satisfies the assigned boundary conditions of free-stream 
velocity at a 3 0° and of zero lift at a”0° and a=90 8 . 
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VELOCITY POTENTIAL * Vx + — tan" ' 


OOWNWASH VELOCITY = - = 

3y 2n x T i + (y/x T ) 2 j 



cje 

— — s 


fcos*^ s«n** + 2 cos 2 * - 4sin ? ot - 3 sin 2 * cos 2 *} 


05 


\ 

! 

\ 


I 


i 


> 


The nonpotential lift equation (A4) and its. >or f vative provide the basis 
for the estimation of lift characteristics for sharp uirfoils and for blunt 
airfoils at transonic and supersonic speeds. The r.i ^tions for zones 3 through 
6 incorporate mult inlying factors to account for • impressibility, sweep, 
aspect-ratio transformations, and detached oow 

Angle of Attack t/p 

The angle of attack for the onset of u-.e leading-edge limit Mach number, 
equation (All), is found by equating equations (A7) and (A10). Equation (A7) 
is the solution to the curvilinear, compressible, potential flow around a 
parabolic leading-edge presented in reference 14. Equation (10) stems from 
the substitution of the limit Mach number (A9) into the general isentropic 
equation (A8). Equation (All) is fitted to various families of airfoils for 
which the leading-edge radius-to-chord ’•atio can be expressed in terms of the 
thickness- to-chord ratio. 



./i + («*7f) 2 (£) 


- 1 . 236 /. / 0.2 + —=• 
V M 2 

_ yi. 528-0.696 M 2 

"( ,v £) 


SET Mg 




(A3, AS) 


(An) 
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Wings of Aspect Ratio Below 2 


The nonpotential lift equations are used for all sharp wings. For aspect 
ratios below 2, the account for aspect ratio differs from the Prandtl trans- 
formation and assumes values taken from reference 15 (fig. 4) and reference 16 
(fig. 6). 


C L - sini*co^«(i-^)^VpiNT)[/AM-(Z-«W)<)tQna 




£ l n 

d0( 


L=-. |)aRWX (?-7$in z o( + 5Sm 4 <x) cos<X+(2 -ARWX)(4 - 10 5irnx+ 5 ^ 


.866 

For Mach numbers above 2, multiply the above equations by ( , 54- ^ - ) 

Equation (A12) for wings alone (ARU*ARUX) can be reduced to the intuitive 
11 leading-edge suction analogy" equation devised by Polhamus (ref, 15), when 
the term — ■ is omitted. 


ry ry 

C- /ARW SINaCOS a+ (2+ARW) SIlTaCOSa 

Jj 


This is the same as equation 4 of reference 15, the sum of the so-called 
potential and vortex lift components. 


C L = Kp Simx co 5^ + K y sin # cos* 


(Ref. 15) 


Camber Drag 


The camber input, CLO, in AEROX is the actual lift coefficient at zero 
angle of attack. The incremental drag coefficient, CDCAM, produced by the 
camber is evaluated on the basis of a constant value of the maximum lift-to- 
drag ratio for the assumed parabolic polar curves for the related symmetrical 
and cambered airfoils. 


CO, 


'CAM “ r~v = 




(A-14) 
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CDC is the minimum drag coefficient of the symmetrical wing, and is one 
of the optional AEROX inputs for each Mach number, when the control integer 
ICDO * 1. When CDO is not specified, and the integer ICDO = 0, CDCAM = CL0/20. 

The parameter CLO, representing wing camber, is attenuated to zero in 
AEROX, when the sonic leading-edge condition is reached. If finite values 
of lift coefficient are to be retained at supersonic speeds, they are input 
as CLOB valuet for each Mach number in the array. (See the shuttle test case 
in Vol. II). 
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APPENDIX B - HORIZONTAL TAIL 
Dovmwash 


The equations describing the dovmwash field used in the TRIMDG routine 
are based on a correlation of theory and data from references 17 through 20. 


Wake angle, EPSD = <X 


3.<I-ARW 


3-AJW 


3/iTtrw 




Nondimensionalized offset distance from tail to wake. 




SPANW 



WhsS. _ JJL\ s/nfeP5D)- sin* (B£) 

5PANW \SPANW ARW / '\ urJW / spanW ^ ) 


Downwash angle at the horizontal tail. 


DWASH = « ^^-ARw'^T 2 /S^ 

[ i/T+t rw 1 / Lt AWJ[ 4^ 


l-i 


2* 


SPANW 


M 


Downwash derivative with respect to angle of attack, d<j/d<x, 


DEPDA* 


3.3-APW 


3-ARW 


^ /It TRW 


«, /^^4^| tcos (<*) 


Supersonic downwash attenuation factor; if 1<M, multiply (B3) and (B4) 
by ATENF , 


ATENF* I 


-o.ipi 


Ton (SWPHILC) 
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Horizontal-Tail Aerodynamics 


> 


The increments in lift coefficient (based on SWING) contributed by an 
aft horizontal tail are evaluated by the nonpotential lift equations. For 
subsonic speeds (Z— 4), 

C L «0R- ^ S.n« „ (l - COS( OW «H)a,5( S 9Hj (BH) 


at supersonic speeds (Z= 5,6), the nonpotential lift for the wing is multiplied 
by suitable ratios of the geometry to obtain the tail lift coefficient. 

C, - C, / ARH \ZaRW-2V COS (SWTVi j SHOR W«,\ rnM , tcin ( B6 ) 

H0R " \ARW/\ARH+2/yC0S(SQH)/^SWINCy\ °^ / V 

01 - <K- DWASH + IT 


The induced drag coefficient contributed by the horizontal tail is positive 
for either up or down tail loads. 



C T TAN(ft+IT) 
HOR 



Pitching-Moment Coefficient and Stability 

The pitching-moment coefficient includes the separate contributions of the 
body, wing and horizontal tail. Initial values of CMO for the wing and CMOJi 
for the body may be input. The wing pitching moment includes the components 
from the upper and lower surfaces and from the chord extension or strake. For 
the transonic limit-shock conditions in flow zone 4, the wing upper-surface 
contributions are further divided into the separate loadings upstream and down- 
stream of the limit shock for both shock geometries covered by IXCD- 0, 1. 
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r r +r k „ . (clwu)(£uMclwl)(el) 
Wr L Mo +LM NOSf^- c&ARW(C 03 fcf 




( c tHo,y^) 

C5ARW(C0Se() 



The offsets EU and EL of the wing surface loadings from the moment center are 
defined for each flow zone Z. 

The static longitudinal stability derivative includes the wing and tail 
components (neglects the body derivative) . 



FLEX and FMEX are the lift and pitching-moment factors for the wing chord 
extension, EBAR is the effective eccentricity of the wing center of lift, and 
LE is the tail length from the reference moment center (usually CG) . 


Longitudinal Trim (ITRIM=1) 


The AEROX program includes an option for evaluating the trimmed values of 
lift and drag coefficients. For angles of attack up to 25°, the trim function 
is performed while maintaining a constant value of lift coefficient. Thus, the 
angle of attack is adjusted to compensate for the changes in tail lift accompanying 
deflection of the tail to trim out the pitching-moment (CM) . 


The horizontal tail deflection to trim while maintaining constant total 
lift coefficient with adjustment of the angle of attack is, 


delh = 


CtA 

f ^m )l lg + cmcl' 

\dOC H J\ C6AXW(C6SK) ) 


UCL1S 
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DCLTR is the increment in lift coefficient produced by the tail deflection DELIi. 
It is recovered through the adjustment in angle of attack, whe^e, 


CX = 
TRIM 


DCLTR 


< ^ + ^^(l-D£PD>rt 

d<X d 0 ( u V ' 


su) 


The incremental drag coefficient DCDTR accompanying the trim process at constant 
lift coefficient is the sum of the increments of induced drag coefficients for 
the wing and tail. 

Cl* -I 


DCDTR* 


(Cl w + dcltr) 2 ^ 

( C L*, R +DCLTR f- C W 

/SwingX 

vARV/ * 

'k'ARH 

\ShoK / 


(i=2) fBIZ) 


D£XP 


DCDTR- - DCLTR (ton oc) + DCLTR [tort + lt+- D£IH)J (3 * c) (B13) 


For angles of attack above 25°, AEROX performs the longitudinal trim at constant 
angle of attack, 0^^= «. 


DCLTR* CM 


( CdARW) 

\ u / 


COSO( 


delh = ~ - r - 


^ ;Jco 5 (dwash) 


DCDTR - 


/dq. 

\d« 

t DCLTR V r * I . 

+ Cos(ohmsn)J " L » 0K ^ £ 
irARij I V 


SWUNG \ 

J* \shor J 


(*•«) 


DCDTR = DCLTR tan(<X-HTf P&lh) 


(36 a) 


Care should be exercised in interpreting the trimmed characteristics. For 
example, if the stability is too large, the tail deflections to trim may exceed 
realistic limits for horizontal tail deflections. The tail angles of attack, 
tfy, should be kept below 40°. 
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APPENDIX C - BODY (BCLCD1 Subroutine) 


The aerodynamic normal force on the body is estimated by the equation 
summing the slender-body contribution and the viscous cross-flow drag. An 
updated discussion of the approach appears in reference 21, which was the 
principal source of the values used for constructing the contour plot of cross- 
flow drag coefficient against cross-flow Reynolds number and cross-flow Mach 
number shown in figure 3. The BCLCD1 subroutine contains explicit equations for 
each of the nine regions indicated- The Reynolds numbers are calculated using 
the atmosphere properties subroutine AT62. 


c Nnose ^[simKcos^fisecT)* cdc(etah) smV jAPLA^^Sm^ 


CLM0S6 - CNNOSf COSH + CLOb 


{Cl) 

'CIA) 


CDNOSe = CNNOSe SM*+ CD OB 


CM N 05 e =• CNMS£ 


/ XC0-0.6XLNU cMnft 
^ CBARW / 


(C1B) 

(CZ) 


At supersonic speeds, X=5,6), the afterbody load is included, for expo 5 ed le*.' ' rh AX a , t . 

CNBOPY * CNM046 (h ( C 3 ) 

CMeoPY - CMN05E XCS) (C4) 


The finite length factor, ETAN, (from a curve fit to ref. 21) is: 


ETAN * 0. 55 +■ 


XLN 

io (BDMAX) 


The slender-body and cross-flow approach formed the basis for the method pre- 
sented in reference 22. 
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